Boundary Lubrication: Squeeze-out Dynamics of a Compressible 2D Liquid 
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The expulsion dynamics of the last liquid monolayer of molecules confined between two surfaces 
has been analyzed by solving the two-dimensional (2D) Navier-Stokes equation for a compressible 
liquid. We find that the squeeze-out is characterized by the parameter go » Po/pc 2 , where Po is the 
average perpendicular (squeezing) pressure, p the liquid (3D) density and c the longitudinal sound 
velocity in the monolayer film. When go <C 1 the result of the earlier incompressible treatment is 
recovered. Numerical results for the squeeze-out time, and for the time-dependence of the radius 
of the squeezed-out region, indicate that compressibility effects may be non-negligible both in time 
and in space. In space, they dominate at the edge of the squeeze-out region. In time, they are 
strongest right at the onset of the squeeze-out process, and just before its completion. 
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I. INTRODUCTION 

Sliding friction is one of the oldest problems in physics, 
and has undoubtedly a huge practical importanc e) 1 ! 2 ' 3 . In 
recent years, the ability to produce durable low-friction 
surfaces and lubricants has become an important factor 
in the miniaturization of moving components in techno- 
logically advanced devices. For such applications, the 
interest is focused on the stability under pressure of thin 
lubricant films, since the complete squeeze-out of the lu- 
bricant from an interface may give rise to cold-welded 
junctions, resulting in high friction and catastrophically 
large wear. 

It has been shown both experimentally and theoreti- 
cally that when simple fluids (quasi-spherical molecules 
and linear hydrocarbons) are confined between atomi- 
cally flat surfaces at microscopic separations, the behav- 
ior of the lubricant is mainly determined by its inter- 
action with the solids that induce layering in the per- 
pendicular directionii^iSiL&iliiSiii. The thinning of the 
lubrication film under applied pressure occurs step-wise, 
by expulsion of individual layers. These layering tran- 
sitions appear to be thermally activatedi 2 * 1 ^*^ 1 ^. Un- 
der strong confinement conditions, some lubricant flu- 
ids become solid-like^SiZ&SiiSiii. Other fluids, notably 
wateipiiiS, remain liquid-like up to the last layer that can 
be removed upon squeezing. This is related to the expan- 
sion of water upon freezingiiiL2&, and should also hold for 
other liquids which expand upon freezing. 

The phenomenology of layering transitions in 2D-solid- 
like boundary lubrication has been studied in ReimS^LS^. 
It has been shown in a series of computer simulations that 
for solid-like layers, the layering transitions are some- 
times initiated by a disordering transition, after which 
the lubricant behaves in a liquid-like manner for the rest 
of the squeeze-out process. Since the typical lateral ex- 
tension in surface force apparatus (SFA) experiments is 



of order 10-100 /zm (much greater than atomic dimen- 
sions) , it is reasonable to expect that during the layering 
transition the squeeze-out can be described in the frame- 
work of 2D continuum fluid mechanics. 

For the first time such squeeze-out layering transi- 
tions were quite recently observed in a chain alcohol, 
C11H23OH (RefiSi 2 ^), by imaging the lateral variation of 
the gap between the two anvil surfaces as a function of 
time. These experiments addressed the n = 1 — > tran- 
sition. More recently, in a refined experimental setup, 
Mugele et al^ were able to image several layering tran- 
sitions (n — ► n — 1, n = 5, 4, 3) of the silicon oil OMCTS 
(spherical molecule, diameter ~ 10 A) in great detail. 
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FIG. 1: Because of the curvature of the solid walls at the 
boundary line, the perpendicular pressure Po will give rise to 
a parallel force component acting on the 2D-lubrication film. 
(Schematic.) 

The basic theory of 2D squeeze-out dynamics was out- 
lined in ReW*. Initially the system is trapped in a 
meta-stable state at the initial film thickness. Squeeze- 
out starts by a thermally activated nucleation process in 
which a density fluctuation forms a small hole, of critical 
radius R c ~ 10 A. Once formed, a 2D pressure differ- 
ence Ap develops between the boundary line separating 
the squeezed out region from the rest of the system, and 
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the outer (roughly circular) boundary line of the contact 
area, thus driving out the rest of the 2D fluid. The origin 
of Ap is the elastic relaxation of the confining solids at 
the boundary line as is illustrated in Fig. ^ 

All earlier analytical studies of squeeze out have as- 
sumed that the lubricant behaves as an incompressible 
2D liquid. While that assumption is quite good for many 
practical situations, recent computer simulations^ have 
shown that, at least at high squeezing pressures, strong 
density fluctuations may occur in the lubrication film. 
For example, in RefiSS it was found that during the lay- 
ering transition n = 2 — > 1, while islands of (temporarily) 
trapped bilayer (n ~ 2) were removed by being squeezed 
into the monolayer, the density of the monolayer film was 
much higher in the region close to the trapped n = 2 is- 
lands than further away. This resulted in a 2D pressure 
gradient in the film which induced a flow of the lubri- 
cant molecules away from the trapped islands. This kind 
of situations clearly calls for a consideration of the finite 
compressibility of the film. 

In this paper we study the dynamics of the expulsion of 
the last liquid monolayer of molecules confined between 
two surfaces by solving the two-dimensional (2D) Navier- 
Stokes equation for a compressible liquid monolayer. We 
find that the squeeze-out is characterized by the param- 
eter go « Pq/ pc 2 , where Po is the average perpendicular 
(squeezing) pressure, p the liquid (3D) density and c the 
longitudinal sound velocity in the monolayer film. When 
go <C 1 the 2D liquid can be considered as incompress- 
ible, in which case the results of the earlier treatment 
are reproduced. We present numerical results for the 
squeeze-out time, and for the time-dependence of the ra- 
dius of the squeezed-out region, for several values of the 
parameter go- The main changes due to compressibility 
occur right at the onset of the squeeze-out process, and 
just before its completion. 



II. THEORY 

We assume the lubricant film to be in a 2D-liquid-likc 
state, and the squeeze-out to be described by the 2D 
Navier-Stokes equations. For a compressible 2D liquid 
these equations take the form 



dn 
~dt 



V • (nv) = 



(1) 



on the nature of the solid walls, e. g. on the amplitude 
of the atomic corrugation and on the structure of the 
solid walls (amorphous vs. crystalline, commensurate vs. 
incommensurate). As long as the lubricant can be con- 
sidered as a 2D fluid, as in the recent squeeze out exper- 
iments by Mugele et al., all the details of the solid walls 
are properly taken into account by the friction coefficient 
fj. In principle fj depends on the normal pressure too, but 
in the measurements of Mugele et al. this dependence is 
slight enough that it appears to be negligible 15 . 

Simple dimensional arguments (see Ref>i&) show that 
one can usually neglect the nonlinear and the viscosity 
terms in (2), and that one can also assume the velocity 
field to change so slowly that the time derivative term 
can be neglected too. Thus, 



Vp + mnrjv = 0. 



(3) 



In what follows we will assume that the squeeze-out 
nucleated in the center (r = 0) of the contact area, 
and spread circularly towards the periphery. The 2D- 
pressure p at the outer boundary r = R of the contact 
area takes the constant value po (the spreading pressure) 
while it takes a higher value p\ at the inner boundary 
towards the n = area (Fig. QJ. In fact, if P(r) is the 
perpendicular pressure acting in the contact area, then 
Pi(r) = po + P(r)a, where a is the thickness of the mono- 
layer (see Ref^). We may in most practical applications 
assume a Hertz contact pressure 



P(r) 



tPo 



r 2 \ 

1_ *0 



1/2 



(4) 



Finally, in order to have a complete set of equations we 
must specify the relation between the 2D pressure p and 
the 2D density n. We will assume that 



p = po + mc 2 (n — no) 



(5) 



where the compressibility is B = 1/mc 2 , c being the lon- 
gitudinal 2D sound velocity. Here no is the 2D lubricant 
density at the periphery r — R of the contact area. 

Writing the velocity field as v = fv(t), equations (1) 
and (3) takes the form 



dn ( d T, 



Q 



(6) 



<9v „ 1 „ 
— + v Vv = Vp 

at ran 



V v + vi VV • v — 77V 



(2) 



where mn and v are the local 2D mass density and the 
velocity of the fluid, p is the 2D-pressure, v\ and are 
viscosities. The last term in (2), i. e. — fjv, describes the 
drag force acting on the fluid as it slides relatively to the 
solid walls' The magnitude of the friction also depends 



dp 
dr 



+ mnnv = 



Using (5) the last equation takes the form 
Combining (6) and (7) gives 



dr 



d 



1\ dn 



dt n \dr ~^ r J dr 



(7) 



(8) 
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The density n satisfies the boundary conditions 

n(R,t)=n (9) 



n(ri(t),t) = n Q 



(10) 



where r = r\(t) is the equation for the inner boundary 
line. Finally, we note that the velocity v(ri(t),t) of the 
2D-liquid at the inner boundary r = r\ (t) must equal the 
radial velocity fi(t) of the boundary line. Thus, if we put 
r = ri{t) and v(ri(t) 7 t) = fi(t) in (7) we get 

— (r 1 (t),t) = -^n(r 1 (t),t)r 1 (11) 

Let us at this point introduce dimensionless variables. If 
we measure the radius r in units of R, time t in units of 
r = fjR 2 / c 2 and density n in units of n Q we get 



dn / d l\ dn 
dt \dr r ) dr 



and the boundary conditions becomes 
n(l,i) = 1 



and 



where 



n{n(t),t) = l+g(t) 



— ( ri (t),t) = -[l + g(t)}h 



9(t) 



-P{n{t)) 



If we assume that P(r) is of the Hertz form, then 



2\l/2 



where 



3o = 



aPo 
mc 2 no 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



Thus, the theory depends only on a single parameter g$. 
Note that no /a « p, where p is the 3D density of the 
liquid, and with the typical values p « 1000 kg/m 3 and 
c ~ 700 m/s we get pc 2 « 500 MPa. In the experimental 
studies by Mugele and Salmeron (as well as in most other 
Surface Force Apparatus studies) the average squeezing 
pressure Pq <C 500 MPa which implies that the liquid 
can be considered as incompressible and the theory de- 
veloped elsewhere can be usecU2*i^i£. However, in many 
practical situations the pressure Pq might be similar to 
the yield stress of the solids which for metals is typically 
of order ~ 1000 MPa. In these cases it is necessary to 



include the finite compressibility of the lubricant in order 
to accurately describe the squeeze-out dynamics. In Sec. 
Ill we shall present numerical results based on the above 
equations for a range of values of the parameter g^. 

Let us first calculate the squeeze-out time in the limit 
when the compressibility B — oo. In this limit, all the 
adsorbates from the region r < r\(t) will be piled-up 
right at the boundary line r = ri . If we consider a small 
angular section A(f> then the driving force acting on the 
boundary line is F — riA(f>[pi —po] = riA(j)P(ri)a. This 
must balance the frictional drag force which is Nmfjri 
where the number of adsorbates N — 7rr^(A0/27r)no. 
Thus we get 



2P(r)a 



(19) 



nomr/ 

Assume first that P = Pq is constant. Thus, (19) gives 

R 2 



2T 



(12) 

where 



T = 



mnQTjR 2 
4P a 



(20) 



(21) 



is the squeeze out time for an incompressible 2D fluid, 
see Refi 12 . Integrating (20) gives 



n(t) 
R 



1/2 



(22) 



so that the squeeze-out time, T*, in the limit B = oo 
is the same as for an incompressible 2D-fluid, T* = T. 
This suggests that the squeeze-out time is independent of 
the compressibility of the 2D-liquid, which our numerical 
simulations presented below indeed show to be the case. 

Next, let us assume that the perpendicular pressure is 
of the Hertz form, appropriate for a fluid between curved 
surfaces. In that case (19) gives 



2T2 



R 



1/2 



or, with r\jR 2 = £, 



/0 (1-C) 1 / 2 ~2T 
Performing the integral gives 



f- 

\R 



1/2 



AT 



Thus, the squeeze-out time for the Hertzian contact pres- 
sure is T]| = AT/3, and the time dependence of ri(t) is 
given by 



n 

R 



St 



2T I V 8T 



4 



ri 
R 



2t 

u 



1/2 



2T* 



1/2 



The squeeze-out time for an incompressible fluid with 
a Hertzian contact pressure is Th = (4T/3)(2 — ln4) « 
0.8183 T which is a factor 2 - ln4 « 0.6137 smaller than 
for the B — oo limiting case. Thus, for a Hertzian con- 
tact pressure the squeeze-out time does depend on the 
compressibility of the 2D lubrication film. 



III. NUMERICAL RESULTS 

We consider first the case of a spatially constant 
squeezing pressure, P(r) = Pq, appropriate for a fluid 
between flat surfaces. In this case numerical calcula- 
tions show that the squeeze-out time is independent of 
the compressibility, i.e., independent of the parameter go- 
However, the time dependence of the squeeze-out radius 
ri(i) does depend on g - In Fig. [Ja) we show this time 
dependence for four cases, go = 0, 0.2, 0.5 and 2. In 
Fig. |2tt>) we show the adsorbate density profile at the 
time point when the squeeze-out radius ri = 0.37?, for 
the same four go -v ahies as in (a). 

Let us now assume the Hertzian squeezing pressure 
profile, P(r) = fk(r). In this case the squeeze-out 
time depends on the compressibility, increasing from 
w 0.8183 T to 1.3333 T as the compressibility increases 
from B — to oo. In Fig. [3] we show the squeeze-out 
time as a function of go- In Fig. 0]we show the time de- 
pendence of ri(t) for four cases, go = 0, 0.2, 0.5 and 2.0. 
In Fig. 0Jb) we show the adsorbate density profile at the 
time point when the squeeze-out radius r\ = 0.37?, for 
the same four go-vahies as in (a). The results of Fig. [21 
and Fig. 0] show that a compressibility parameter go as 
small as 0.2 should produce a measurable difference in 
the squeeze-out evolution at all times. The main effect 
of compressibility appears at the beginning and at the 
end of the squeeze-out process. Initially, compressibil- 
ity favors piling up of fluid at the squeeze-out bound- 
ary, which can as a result expand more rapidly, com- 
pared with the case of an incompressible fluid. On the 
other hand, when the hole approaches the boundary of 
the contact region, the squeezing-out speed of the com- 
pressed fluid is smaller, due to its increased density and, 
consequently, friction. In the case of uniform squeezing 
pressure these two effects compensate exactly, leading to 
a total squeeze-out time independent of the compressibil- 
ity. For a Hertzian pressure distribution, the squeeze-out 
time instead increases with increasing compressibility: 
the initial speed up is overcompensated by the enhanced 
friction at the periphery of the contact area. 
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FIG. 2: Constant squeezing pressure, P(r) — Po- (a) 
Squeeze-out radius r\ (in units of the radius R of the contact 
area) versus time (in units of the squeeze out time). We show 
results for four different cases, namely go = (corresponding 
to an incompressible adsorbate layer), 0.2, 0.5 and 2.0. (b) 
Adsorbate density distribution n(r) (in units of the natural 
density no) during squeeze-out, for the same four cases as in 
(a) at a time when the squeeze out radius n = 0.37?. 




FIG. 3: Dependence of the squeeze-out time on the com- 
pressibility parameter go for a Hertzian squeezing pressure. 
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IV. SUMMARY 



The continuum mechanics theory of squeeze-out has 
been solved numerically for a 2D compressible liquid. We 
considered both a constant normal stress, and a Hertzian 
normal stress, and assumed a centro-symmctric squeeze- 
out. We found that the squeeze-out is completely char- 
acterized by the compressibility parameter g w P /pc 2 , 
where Pq is the average perpendicular squeezing pressure, 
p the (3D) liquid density and c the longitudinal sound ve- 
locity in the monolayer film. When go <C 1 the 2D liquid 
can be considered incompressible, and the earlier results 
are reproduced. We presented numerical results for the 
squeeze-out time, and for the time dependence of the 
squeeze-out radius, for a grid of values of the parameter 
go- For a constant squeezing pressure, the squeeze-out 
time was found to be independent of the compressibility 
parameter go, while for a Hertzian contact pressure it in- 
creased slightly with increasing go. It is hoped that these 
theoretical results will soon be submitted to experimental 
check. 
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FIG. 4: Hertzian squeezing pressure, P(r) given by Eq. (4). 
(a) Squeeze-out radius n (in units of the radius R of the con- 
tact area) versus time (in units of the squeeze out time). We 
show results for four different cases, namely go = (corre- 
sponding to an incompressible adsorbate layer), 0.2, 0.5 and 
2.0. (b) Adsorbate density distribution n{r) (in units of the 
natural density no) during squeeze-out, for the same four cases 
as in (a) at a time when the squeeze out radius n = 0.3R. 
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